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We consider a Quantum Computer with n quantum-bits ('qubits'), where each qubit is coupled 
independently to an environment affecting the state in a dephasing or depolarizing way. For mixed 
states we suggest a quantification for the property of showing quantum uncertainty on the macro- 
scopic level. We illustrate in which sense a large parameter can be seen as an indicator for large 
entanglement and give hypersurfaces enclosing the set of separable states. Using methods of the 
classical theory of maximum likelihood estimation we prove that this parameter is decreasing with 
l/\/n for all those states which have been exposed to the environment. 

Furthermore we consider a Quantum Computer with perfect 1-qubit gates and 2-qubit gates with 
depolarizing error and show that any state which can be obtained from a separable initial state lies 
inbetween a family of pairs of certain hypersurfaces parallel to those enclosing the separable ones. 



I. INTRODUCTION 

The sensitivity of quantum systems to interactions 
with the environment is one of the challenging prob- 
lems for the realization of Quantum Computers j|lp||. 
But apart from this motivation the effect of the envi- 
ronment to quantum states has been subject of pure re- 
search for many decades (see Q and references therein). 
After all the decoherence caused by the environment is 
commonly accepted to be the explanation for classical 
behavior of physical systems in every-day life i.e., 
on the macroscopic level of physics [QJ^] . Although there 
is no precise definition of the word 'macroscopic', most of 
those explanations contain (explicitly or implicitly) the 
statement that the destruction of coherence takes place 
on a very small time scale particularly for superpositions 
of 'macroscopically distinct' states, i.e., States showing 
quantum uncertainty on the macroscopic level. Despite 
the fact, that this statement cannot be maintained with- 
out taking into account the way of coupling to the envi- 
ronment (see Q), it has served as an intuitive motivation 
for our investigation of the sensitivity of many particle 
quantum states with respect to a coupling to independent 
environments. For this we introduce a function e : p h> e p 
from the set of states to the positive numbers quanti- 
fying the property of showing quantum uncertainty on 
the macroscopic level and prove quantitative statements 
about the sensitivity of those states p having large values 

We show that large values for e p require large-scale- 
entanglement in the sense that there cannot be small 
clusters of entangled qubits without entanglement be- 
tween qubits in different clusters. Therefore our in- 
vestigations should be considered in the context of the 
fragility of entanglement, which is an important subject 
since it is decisive for the computational power of deco- 



hered Quantum Computers [l0| . 

The connection between the sensitivity of a state with 
respect to disturbances of the environment and the prop- 
erty of being a superposition of macroscopically distinct 
states can be illustrated by the following straightforward 
example: Take an n-qubit Quantum Computer, that is a 
quantum system with the Hilbert space 

C 2 ® . . . ® C 2 . 
, ' 

n 

Denote the canonical basis states by the binary words of 
length n. Furthermore we use the following terminology: 
The Hamming weight of a binary word is the number of 
characters T'. The Hamming distance of two words is 
the Hamming weight of their difference. Now we take a 
superposition of two arbitrary basis states 

m :=^(|a> + |&». 

If a and b differ at many positions, i.e., have a large 
Hamming distance, we call \a) and \b) macroscopically 
distinct states. Now we perform a measurement of in 
the canonical basis |0), |f) of C 2 on one randomly chosen 
qubit i. Obviously this state collapses to \a) or |6) if and 
only if the words a and b differ at the position i. Hence 
the superposition state is more fragile if a and b have a 
large Hamming distance. At the first sight it seems to 
be straightforward to characterize the fragility of a state 
(with respect to dephasing) by the probability that it is 
changed by a measurement of a randomly chosen qubit. 
But this probability is f in the generic case: Even the 
unentangled state 

k) := (i=(|0) + |l»)®" 
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is changed by a measurement of any qubit 1 . Nevertheless 
we want to consider this state as much less fragile than 
the 'cat-state' 

It) :=^(|o...o) + |i...i» 

since the 'error' caused by the single-qubit-measurement 
can be corrected by a single-qubit-operation in the case of 
the state \tt) whereas the local disturbance of the cat state 
1 7) requires a much more complicated procedure restor- 
ing the entanglement. Therefore we consider fragility as 
a property of a class of states rather than of a single 
state. We prove a class of highly entangled states to be 
fragile in the sense that every mixed state obtained by 
small independent disturbances of each qubit lies outside 
this set. 



II. TWO SUFFICIENT CONDITIONS FOR 
ENTANGLEMENT IN MANY PARTICLE 
SYSTEMS 

As usual we call a state of the Quantum Computer 
a product state if its density matrix is an n-fold tensor 
product of the form 

Pi <E> P2 <E> • ■ • <8> p n 

where each pi is the density matrix of the qubit i. A sep- 
arable state is a convex combination of arbitrary many 
product states. For any operator a acting on an Hilbcrt 
space we denote its operator norm by ||a||. By an '1- 
qubit-operator at the qubit V (or 'acting on the qubit V) 
we mean an operator acting on the n-fold tensor product 
of C 2 which is of the form 



l®l®...®a<g>. 



.1(g) 1 
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for any a / 1 acting on C . 

Furthermore we introduce the following type of ob- 
servables: Let (<2j)i<i<„ with ||a.j|| < 1 be a family of 
self adjoint operators where each a, is acting on the qubit 
i. Then we define the averaging observable 



ri — ' 



In the case that the difference of the lowest and the 
greatest eigenvalue is the same for every a% we call a 
an equally weighted average of 1-qubit observables. De- 
spite the fact, that there is no precise distinction between 



the macroscopic and microscopic level in an n-qubit sys- 
tem, we have good reasons for considering the equally 
weighted averages as the 'most macroscopic ones'. An 
easy example might illustrate this point of view: If the 
qubits are represented as spin-l/2-particles, the mean- 
magnetization of the system in z-direction is given by 
the averaging observable obtained by setting 



■= a® 



where is the copy of the Pauli matrix a z acting on 
the qubit i. If we would define 

(H := AjO-W 

with arbitrary Aj we get a less macroscopic observable in 
general since a is dominated by the spins of those i with 
large |A 4 |. 

In a product state there is no correlation between 
the values of two 1-qubit observables at different qubits. 
Hence we conclude the following from classical probabil- 
ity theory: 

For any observable a and any density matrix v let s av 
be the standard deviation of a in the state v, i.e., 



s a , v ■= \Jtr{ya 2 ) - {tr(ya)) 2 . 

Let p be a product state. Then for the standard devi- 
ation of any averaging observable a in the state p the 
inequality 
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(1) 



holds since the variance of a sum of independent random 
variables is the sum of their variances and the variance 
of any at cannot exceed 1 due to its operator norm. We 
conclude, that every separable state has a decomposition 
into states fulfilling inequality (Q) for every family (aj). If 
we want to use this result for showing that a given state 
is not separable one would have to check every possible 
decomposition into pure states. Hence one might doubt 
its practical importance. However, we can derive another 
sufficient condition for entanglement which does not re- 
quire to check every decomposition: The question, as to 
which extent a state can be decomposed into those pure 
states with small standard deviations with respect to any 
given observable a, is closely related to the question as 
to which extent its density matrix (written in any basis 
diagonalizing a) is dominated by 'strongly' off-diagonal 
positions, i.e., those positions where row and column cor- 
respond to rather different eigenvalues of a. The convex 
function 



x In this state is taken as an example for the difficulty of 
maintaining coherence in large Quantum Computers. 
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p i > sup \tr(p[a,b])\ 

\\b\\<l 

can be considered as a measure for the 'dominance of the 
strongly off-diagonal' terms since it vanishes for every p 
commuting with a. 

To be precise, we have the following lemma: 

Lemma 1 Let A be an arbitrary (finite dimensional) 
matrix algebra. Let b,c £ A with b selfadjoint and 
\\c\\ < 1. Let v G A be an arbitrary density matrix. Then 
we have the following inequality: 

\tr(v[b,c])\ < 2s blV (2) 

Proof: Since the standard deviation is a concave func- 
tion on the set of probability measures on R, we have 
Sb.w > J2j ^j s b,uj if v is the convex sum v := J2j ^j v j- 
Therefore we can assume v to be a density matrix of a 
pure state, i.e. v — |^>)(i/>|. Let us choose an eigenvector 
basis of b and expand c and \t\)) with respect to this ba- 
sis. Let Ai, .., Afc be the set of eigenvalues of b with their 
corresponding multiplicities and ipj be the coordinates of 
\ip). We have: 

M[b,c]i>)\ = \J2M*i-^)^ j \ 

< | V'iCAi - M)cii^il + I (Aj - n)ipj\, (3) 

ij ij 

for any /i£R. 

Defining the vector n with Ki := (Xi — p)ijji. the in- 
equality (Q) reads 

|<V#,c#)| < \(k\oI>)\ + \{iI>\ck)\, 

which gives us an upper bound by the Cauchy Schwartz 
inequality 

!<#&,<#> I < 2||«|| ||c</>|| <2||«||, 

where the last estimation holds due to the operator 
norm of c. With the definition p, :— (ip\bil>) the vector 
norm ||k|| is the standard deviation of the observable b. 
□ 

By the triangle inequality we conclude: 

Corollary 1 Let p be a density matrix of a finite dimen- 
sional quantum system. Let b be an arbitrary selfadjoint 
operator and c be an arbitrary operator with \\c\\ < 1. 
Assume p to have a decomposition into pure states of the 
form 

P = Y1 X JP^ ( 4 ) 

3 

where Xj > and • Xj = 1 . The states pj are arbi- 
trary pure states. Then we call a := Sb, Pj 'the mean 
standard deviation of the observable b in the state p with 



respect to the decomposition (yj ' and have the following 
inequality: 

^\tr{p[b,c])\<a 

Note that for any pair of selfadjoint operators b, c and 
any state v the expectation value i tr(v[b, c]) is real. For 
any such pair b, c we define the hypersurface 

H b ,c,r ■= {p G S | i tr(p[b, c]) = r} 

where S is the set of density matrices of the n-qubit sys- 
tem. We conclude from Corollary [t| 

Corollary 2 ('Hypersurface- Criterion') Let a be an av- 
eraging observable as in the beginning of this section and 
c be a selfadjoint operator with \\c\\ < 1. Then every sep- 
arable state p lies inbetween the hypersurfaces H- C:i __2_, 

i.e., 

This gives us a sufficient condition for entanglement 
which is easy to handle since it can be verified by finding 
just one pair a, c such that inequality || is violated. 

Despite the fact, that there is no commonly accepted 
quantification of entanglement (see |ll]-[l3|] ) , we will con- 
sider a large value (compared to 1/ y/n) of the term 
tr{p\a,b]) for any such family (a^) and any such b as a 
sufficient condition for p to be 'highly entangled'. The 
term highly entangled might be interpreted in two differ- 
ent ways: Firstly a large value shows that the state has 
a great distance from the set of separable states in the 
trace norm. Secondly it shows that there is entanglement 
between many qubits: 

Lemma 2 Take a partition of the qubits {1, . . . , n} into 
subsets ('clusters') of size l\, . . .,lk with the property that 
the state p has no entanglement between qubits of differ- 
ent clusters. Let a be an averaging observable and b be 
a selfadjoint operator with \\b\\ < 1. Then the following 
inequality holds: 

\tr(p[a,b])\<l (6) 

Y i<k 

Proof: Denote the clusters by Si , . . . , Sk C 
{1, . . . , n}. It is sufficient to take a state which is fac- 
toring with respect to this partition into clusters since 
the set of states fulfilling inequality (^) is convex. For 
such a 'partial product state' p the standard deviation of 
a is given by 
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where s, denotes the standard deviation of the observable 
Sjes a i wn i cn is l ess or equal to i,-. Lemma [l] completes 
the proof. □ 

In the following, we will restrict ourselves to the 
equally weighted averaging observables. Up to a constant 
factor and a constant summand they can be obtained by 
taking every <Zj as a projection. 

Therefore, we shall consider the parameter 

e p :=sup|tr(H£,&])| (?) 

Q.b 

where every Q is the average Q over the 1-qubit pro- 
jections Qk and ||6|| < 1, as a reasonable quantification 
for the property of showing quantum uncertainty on the 
macroscopic level. 

Example 1 In order to give more intuition about the 
states with large e p , we assume p to be the density matrix 
of a coherent superposition of two distinct basis states \f) 
and \g), i.e., p = ±(|/) + \g))((f \ + (g\). Let f and g be 
binary words with Hamming weights wgt(f) and wgt(g). 
Let Pk be the projection onto the state |1) for the k-th 
qubit. With the definition 

b:=i\f){g\-i\g){f\. 

we get 

tr(p[P,b}) = l -{wgt(f) - wgt(g)), 

and hence we have e p > — \wgt(f) — wgt(g)\. 

We prove the following more general statement: 

Lemma 3 Let fi , . . . fj , g\ , . . . , Qj be a set of 2j distinct 
binary words. Let p be the density matrix given by 

p := ^2\k\ipk)(ipk\ 
k 

where X k is the probability of the pure superposition state 
m := ±={\f h ) + \g k )). 

Define b by 

b:=Y.^)^\-^\g k ){h\). 

k 

Then we have the following equation: 

tr(p[P, b}) =i~y2 ^k(wgt(f k ) - wgt(g k )). 
k 



Proof: Using 

P\f) = -wgt(f)\f) 

n 

for every binary word /, the statement follows by easy 
calculations. □ 

Note that the operator b in the definition above fulfills 
the requirement ||6|| = 1 since the operators 

i\fk)(gk\ ~ i\gk)(fk\ 

have operator norm 1 and act on mutually orthogonal 
subspaces. 

The states \ip k ) are superpositions of macroscopic dis- 
tinct states if the difference wgt(f k ) — wgt(g k ) has the 
order of n rather than the order of 1. In this case we say 
the state shows quantum uncertainty on the macroscopic 
level. The sum *^2 k X k {wgt{f k ) — wgt(g k )) measures to 
what extend the mixture p consists of pure states with a 
large uncertainty of the observable P. Note that the mix- 
ture of two states with large parameter e can have small 
e due to the fact that the mixture of two highly entan- 
gled states can be separable. Therefore the assumption 
that the 2j binary words f\ , . . . , fj , gi , . . . , gj are mutu- 
ally distinct is essential and turns up not to be just a 
technical requirement for the proof: Take the mixture 
given by 

P'-= ^(1^1)^x1 + IVaXV*!) 

with 

|W 2 ) :=^(|0...0)±|l...l)). 

Easy calculation in the canonical basis shows that 
tr(p[P, b]) vanishes for every operator b since p and P 
are diagonal in this basis. Actually, we can get this re- 
sult using Lemma |l| as well: The state p has another 
decomposition into the pure states po := |0 . . . 0) (0 . . . 0| 
and pi :— |1 . . . 1) (1 . . . 1|. Both states do not show any 
uncertainty with respect to the observable P, i.e., the 
standard deviations s-p „ and s-p „ vanish. 

We can make general statements about the range of the 
convex function p i— > e p on the set of density matrices: 

Lemma 4 For every n the range of the function p t— > e p 
is the interval [0, 1]. 

Proof: If p is the maximally mixed state, i.e., p is the 
identity matrix up to a constant factor, we have e p = 
since tr(p[a, b)) = tr{[a, b}) = for every pair of operators 
a, b. In general we have e p < 1 due to Lemma [j] since 
the standard deviation of P cannot exceed 1/2 due to 
the fact that its spectrum is contained in [0, 1]. For the 
cat state as defined in the introduction we can conclude 
£p > 1 by setting / := . . . and g := 1 . . . 1 in Example 
pi Hence for the cat state we have e p = 1. We can obtain 
any value between and 1 by a mixture of the cat state 
and the maximally mixed state with the corresponding 
weight. □ 
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III. THE ERROR MODELS 

In many cases it is well-justified from a physical point 
of view to assume errors acting independently on every 
qubit (l(J . One kind of these 1-qubit-error which seems 
reasonable is a random dephasing with respect to the 
canonical basis. Describing this on the set of density 
matrices this error affects the state like a measurement 
instrument. We describe the effect of the dephazing en- 
vironment by a map G on 2 the set of n-qubit density 
matrices as follows: 

Let Pi be as in Example ^ and Mi the instrument per- 
forming a measurement of the qubit i in the canonical 
basis of C 2 , i.e., 

M i (j>) = P i pPi + (l-P i )p(l-Pi) 

for every density matrix p. Then our first error model 
will be an instrument G which acts on each qubit inde- 
pendently as wMi + (1 — w)id, i.e. G acts as: 

G := JJ(u;Mi + (1 - w)id) 

i 

where id denotes the identity map and w is the error 
probability. 

Our second error model is given by depolarizing chan- 
nels acting on each qubit independently: Let I be this 
map on the set of 1-qubit density matrices which maps 
every state to the maximally mixed one. Let Ii be the 
canonical extension of this map from the state space of 
the qubit i to the n-qubit-system, i.e., 

Ii := id ® . . . ® id ® I ® . . . ® id. 

Define the instrument D by: 

D := Y[(wli + (1 - w)id), 

i 

In the following chapter we shall study the images of the 
maps D and G according to the error probability w and 
the size n. 



( ^ j LeCi i£L 

Then easy calculation shows that G can be written as 
the convex combination of all the Gi with binomial coef- 
ficients: 

n 

G = Y J B nw (l)G l . (8) 

where we use the abbreviation 

B nw {l) := (j^j w l (l~w) n - 1 . 

The instrument Gi is a random machine performing a 
measurement on every qubit in a randomly selected 1- 
element subset of qubits. It can map a pure state to a 
mixed one for two reasons: Firstly we do not know, which 
/ qubits are measured, i.e., which L was selected, and sec- 
ondly we do not know the measured result. If we knew 
both, we would get a certain pure state obtained from the 
original one by a partial collapse of the wavefunction. In 
the following we show, that most likely this collapse leads 
to a state in which the observable P = - V" P; has small 
standard deviation provided that I » 1. Intuitively this 
is not astonishing since the measurement of / qubits al- 
lows a prediction of the values of the average observable 
P with a high 'confidence level'. This analogy to the the- 
ory of maximum likelihood estimation motivates the idea 
of the proof in the following quantitative analysis: 

Theorem 1 Let p be an arbitrary density matrix of an 
n qubit system. Let Gi and P as described above. For b 
an arbitrary operator on (C 2 )®" with \\b\\ < 1 we have 



IV. QUANTITATIVE STATEMENTS ABOUT 
FRAGILITY 

In order to investigate the way in which the instrument 
G affects a state we introduce a family of instruments 
{Gi)i< n which is defined as follows: Let £; be the set of 
1-element subset of {1, . . . , n}. Set 



Before we prove the theorem we draw some conclu- 
sions. In order to get statements about G(p) instead of 
Gi (p) we can use the decomposition (0) . For this we have 
to consider the case I = separately. The instrument Go 
is the identity map and occurs in the sum ^) with the 
weight (1 — w) n . The standard deviation of P can never 
exceed 1/2. Therefore, \tr{G a {p)\P ,b])\ = \tr(p[P,b})\ < 
1 and we conclude: 



2 G is a completely positive trace preserving map (see JlJ]), 
since every manipulation of a quantum state can be described 
by a map of this type. 
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Corollary 3 Let G be as in equation (j<J). Then we have 
for any arbitrary density matrix p 



\tr(G(p)[P,b})\ • Z~T^ !< " 



v^Vn-l' 



+ (i-wy 



Proof (of the Theorem): Since the set of density 
matrices fulfilling the inequality is convex, we may re- 
strict the proof to the case of p being a pure state, i.e. 
p = with \ip) £ (C 2 )® 71 . Let us assume a mea- 

surement which has been performed on the qubits in L, 
where L is an arbitrary /-element subset of {1, . . . ,n}. 
Let Rl be the map Rl ■ {0, 1}™ — * {0, 1}' restricting a 
binary word to the set L. Let Ph,g be the projector onto 
the linear span of those basis vectors given by the binary 
words in R~^ l {g) for any g 6 {0, 1}'. 

Then \tp)(ip\ is transduced to the mixed state 

£ PL, g ma>\p L , g . 

gG{0,l}> 

On the set Ci we introduce the measure r/ as the 
equally distributed probability measure, i.e. 



VLG C, : n(L) 



1 



Then G; transduces \ip)(ip\ to the state 

Gi(W(4>\) = r '( x ) E p LMQ>\PL, g . (io) 

LeCt ge{osy 
With the definition \ij) L<g ) := \P L , g i') /\\PL, g ip\\ and 



p(L,g) := ri (L)\\P L , g ^ 



(11) 



we obtain: 



Le£ !:S e{o,i} ! 

Note that p(L,g) is the probability for the event 'a 
measurement has been performed on the qubits in L and 
the result (<?i,.. . ,gi) (in an ascending order) has been 
obtained.' We denote this event by (L,g). 

Let SL, g be the standard deviation of P in the state 

\ll>L,g){lpL,g\- 

In order to prove the theorem it is sufficient (see Corol- 
lary [I]) to show 



E^, 5 K, 9 <^V^-T 



(12) 



Ci x {0, 1}" 

endowed with the product measure r; ® <jy, , where as- 
signs to every binary word the probability given by the 
square of the probability amplitudes of %p. All the random 
variables which will be introduced below are defined on 
this product space. In the formal framework of this space 
the formal correct notation of (L, g) is {L} x i?^ 1 (g). We 
will keep the less formal notation (L,g) for reasons of 
convenience. Now we introduce the random variable 



by 



U : Ci x {0,1}" -► R 



U(L,b) := -wgtfb). 
n 



Then we can write the standard deviation in inequality 
© as 



sl,, 



y/EttU-E(U\L,g))*\L,g), 



where E(.\L,g) denotes the expectation value of a ran- 
dom variable with the conditional probability measure 
given the event (L,g). Furthermore we define the ran- 
dom variable S by 

S(L,b) := jwgt(R L (bj). 

Due to the fact, that S has a constant value on every 
subset (L,g), we can give the following upper bound: 



s L , g < y/EttU-S)*\L,g), 

because for any arbitrary random variable X and /i£E 
the expectation value E((X — p) 2 ) is minimized by p — 
E(X). 

Since the square root function is concave we get 



^p(L >9 )s L , g < Y,P^9)E((U-Sy\L,g) (13) 



= / E P(b)E((U - Sr\b) (14) 
Y be{o,i}" 

The last equation holds since the family of sets 

(L,g)Lec ug e{a,iy 

as well as b e {0, 1}™ define different partitions of the 
probability space & x {0, 1}". We conclude 



For doing so we introduce the probability space 



Vp(i,«)si, s < sup y/E((U - S) 2 \b). (15) 

L,g b ^ ^ 

Note that for any fixed b the product SI is a random 
variable with a hypergeometric distribution, since it mea- 
sures the hamming weight of the restriction of b to a 



G 



randomly chosen Z-element subset of {l,...,n}. Fur- 
thermore U is a constant for any given b and U is the 
expectation value of the random variable S with respect 
to the conditional probability measure given the event 
b. Therefore the square root in the right hand term is 
the standard deviation a of S and is given by ( jlq] . Sec. 
2.3): 



I 



l{n-l) 



pq 



with p :— wgt(b)/n and q := 1 — p. Since pq < 1/4 
independent of wgt(b) we can estimate the term ( |l5| ) by 



n - 1 



□ 



Since G and G; describe error models which are not 
invariant with respect to local unitary transformations, 
there is no evident generalization for any other family 
(Qi) of projections instead of (Pi). In contrast, the error 
model defined by the map D (see end of section i s 
symmetric with respect to the group SU2®SU2>S)- ■ .&SU2 
of independent local unitary transformations on every 
qubit. Therefore we obtain estimations for the standard 
deviations of any observable obtained by averaging over 
an arbitrary family (Qi) of 1-qubit-projections and get: 

Theorem 2 Let p be an arbitrary state of an n-qubit 
Quantum Computer. Let r wn as in (j^). Then we have 
the following inequality: 

e D(p) < T wn . 



Proof: Let a x ^ be the operator representing the Pauli 
matrix o~ x acting on the qubit i. Let M; be as in sec- 
III and Fi be the instrument performing the bit-flip 



tion 



v i — > o~ x ' vo~x ' on every density matrix v. 
Then 7, of section [II is given by: 



h = (^(Fi + id) o Mi). 
Hence D can be written as the product 



D 



IlA, 



where Di is defined by: 

D % = i(^( F i + id ) ° M * W ) + i 1 ~ w ) id )- 

In analogy to equation (||) we can decompose D into a 
convex sum of instruments D\ where Di is a machine 
performing a depolarizing error on a randomly chosen l- 
element subset L C {1, . . . , ra} of qubits. Hence we have: 



A(IWI)= E r,(i)nJ(^ + «)^iW)M) (16) 
Led ieL 

= E r l (L)Hl(F l + l d)Y[MmW)- (17) 



L£C t 



i£L 



i£L 



Using 



l[M i (\i>)M)= E \\Pl,M 2 W>l,9)^l, 9 

g£{0,l}' 



i£L 



and the definition (fL l|) we obtain 



e p(^fl)n^ F <+ id )^o<^, ff i) as) 

LeCi ,g£{0,l} 1 i( ^ L 

e p(i,5)Ein^(i^)(^D ( i9 ) 



LGCi,g£{0,l}> 



2 l 

T<ZL i£T 



This completes the proof: For any T the standard de- 
viation of P in the pure state Yl ieT Fi(\ipL,g)(i>L, g \) is 
the same as in the state \ipL, g ){ipL, g \, since every Fi is 
only a permutation of the states |0) and |1) in the qubit 
i. 

Therefore we have shown, that Di(\ip)(ip\) has a de- 
composition such that the corresponding mean standard 
deviation of P is less or equal to 



1 



- I 



Z(n-l) 



Therefore the mean standard deviation of D(\ifj)(ip\) is 
less or equal to r wn /2 (see the Definition of r wn in (||) 
and the decomposition of D into a convex sum of D{) 
and hence 



\tr(D(p)[P,, 



< r v 



by Corollary 111. Due to the symmetry of the error 
model with respect to local unitary transformations we 
can substitute P by any other average Q over 1-qubit- 
projections. □ 

The theorem shows, that the 'entanglement- 
parameter' e p of any state p is extremely sensitive to 
small depolarizing perturbations acting on every qubit 
independently. 

At the first sight, the only quintessence of these results 
seems to be that they support the well-known fragility 
of entanglement by a quantitative analysis without tak- 
ing into account the possibility of error correction [ fl6| . 
However, we can make easy conclusions for the following 
model which is so general as to include every possible 
error correction procedure. We take a Quantum Com- 
puter which allows only an imperfect implementation of 
gates. Taking into account that every error correction 
has to rely on these gates producing new errors, we see 
that there arc states which never can be obtained. 



7 



The following theorem shows this statement quantita- 
tively. Therefor we assume that our Quantum Computer 
is endowed with the following set of operations on the set 
of density matrices: 

• perfect 1-qubit gates, i.e., maps of the type p i— > 
upu* where u is an arbitrary unitary operator act- 
ing on 1 qubit and p is the density matrix of the 
Quantum Computer 

• imperfect 2-qubit gates g of the following type: 

9{p) ■= (1 - w)upu* + w(Ii o Ij)(p), 

where w is the error probability and u is an arbi- 
trary unitary operator acting on the qubits i and 
3- 



Theorem 3 Let the Quantum Computer ( endowed with 
the basic operations above) be initialized in a separable 
state. Let (Qi) be a family of projections where Qi is 
acting on the qubit i. Let b with ||6|| < 1 be a self adjoint 
operator and r wn as in Corollary y. 

Then it is not possible to prepare a state outside the 
slice described by the pair of hypersurfaces 



b,±x 



with 



x := — sup{r wk k + \Jn — k}, 

n k<n 

i.e., it is not possible to prepare a state p with e p > x. 

Proof: Let the initial state be a pure product state. 
Let k be the number of those qubits which are accessed by 
2-qubit-gates during the preparation procedure. Without 
loss of generality assume them to be the qubits 1, . . . , k. 
Let p be the state obtained by the preparation. Then we 
have: 

k 1 1 ™ 

\tr(p[P,b])\<-\tr(p[-J2^b])\ + -Mp[ ]T Pi,b})\- 

n k * — » n t — * 



i<k 



i=k+l 



Firstly we can show 



(20) 



i<k 



by a slight modification of the arguments in the proof of 
Theorem |l|: Due to the fact, that for any qubit i with 
i G {1, . . . , k} there is a step in the algorithm which is the 
last access on i by 2-qubit-gates, every of those qubits is 
exposed to the depolarizing channel, i.e., the state p can 
be written as p = D(v) <B> p., where v and p are states of 
the qubits 1, . . . , k and k + 1, . . . , n, respectively and D 
is the fc-fold depolarizing channel on the qubits 1, . . . , k. 
The map D is given by the restriction of Y\ i<k Di to this 



subsystem. By convexity arguments, it is sufficient to 
show 



IM((^)®^)iE p - 6 ])l 



< r„ 



i<k 



for every pure state v. Since D is the analogue to D 
for the fc-qubit-system 1, . . . , fc, the mean standard de- 
viations of the mixture Dv can be estimated in analogy 
to the proof of Theorem [l]. Using Corollary [j] shows in- 
equality (po|). 

Since the qubits k + 1, . . . ,n are still in a product 
state after the preparation procedure, the observables 
{Pi]i>k+i are stochastically independent. Therefore the 
variance of the sum J2i>k+i ^» ^ s * ne sum 01 * ne variances 
of Pi. Hence the standard deviation of ^ X)i>fc+i Pi can " 
not exceed — k. Hence we have 

\tr(p[- y P l: b})\ < -Vn~k. 

i>k+l 



The extension of the proof from a pure to a separable 
initial state is given by standard convexity arguments. □ 

Corollary 4 Modify the assumptions of Theorem as 
follows: Instead of 2-qubit-gates with depolarizing error 
we assume a dephasing error, that is, we have imperfect 
2-qubit-gates g of the form 

g(p) = (1 - w)upu* + w(Mi o Mj)(p), 

where u is an arbitrary unitary acting on the qubits i and 
3- 

Then it is not possible to prepare a state outside the 
slice defined by the hypersurfaces 

TT 

P,b,±x 

with b and x as in Theorem [|. 

Proof: In analogy to the proof of Theorem |[ Those 
qubits 1, . . . , fc, which are accessed at all, are subjected 
to the dephasing error, i.e., to the instrument G on the 
fc-qubit Quantum Computer 1, . . . , fc. □ 

Theorem |^ shows up to which values the error rate 
w of a depolarizing channel has to be reduced in order 
to maintain macroscopic superpositions in the sense dis- 
cussed here: In order to investigate the asymptotic be- 
haviour of 



n 1 



1=0 



Vl V 71-1 



B nw (l) + (l-w) n 



we choose an arbitrary < a < 1 and split the summa- 
tion over I into 2 parts and get 
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El In — I „ 



VI V n - 1 

KKwan 



(21) 



+ E 7f\/|^T B ™ (0 + (1 ~ wri (22) 

Q!'i/jn<Z<n 

< V S„„(0 + -^= + (l-w) n (23) 



KKwna 



— ^ B nw (l) H == 



Q<l<wna 



s/nwa 



Due to the Tschebyscheff inequality the sum over the 
binomial coefficients from to wna is less or equal to 
1 Hence we conclude for any < a < 1: 



l(l-Q) a 



1 



l 



ii>7i(l — a) 2 ^nwa 



V. CONCLUSIONS 

We have shown, in which way the parameter e p and the 
'hypersurface-criterion' can be used for detecting many- 
particle-entanglement. States with large e p are fragile in 



the following sense: We see that r wn is decreasing with 
0(y/wn) for n going to infinity. For a fixed error proba- 
bility w we conclude that the maximum of e p which can 
be attained by states subjected to the n-fold depolarizing 
channel is decreasing with 0(l/y/n). 

Furthermore we see that the preparation of states p 
with a given e p by imperfect gates requires an error prob- 
ability decreasing with 1/n or faster. Similarly, it re- 
quires an error probability not greater than 0(1/ n) in 
(24) order to have a state with a fixed e p after one time step 
if the decohering effect of the environment during one 
time step is described by the map D. Note that these 
are statements about the physical states of the Quantum 
Computer in contrast to the logical states which may be 
defined by a certain Quantum Code. 

At the moment, we cannot see whether stronger 
bounds can be given. However, it should be empha- 
sized that our bounds for the fragility of macroscopic 
superpositions are much weaker than those which are 
suggested by simple but non-generic examples: The in- 
strument D transduces the density matrix of the cat state 
^75 (|0 . . . 0) + 1 1 ... 1) ) to a density matrix p containing 
the cat state with probability (1 — w) n and containing un- 
cntangled states with probability 1 — (1 — w) n . Therefore 
we get exponential decay of e p for increasing n. 
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